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•  Extended  earlier  work  on  aluminum-water  burning  underwater  propulsion  system. 

•  Mass  estimation  methodology  extends  energy  density  estimates  for  neutral  buoyancy. 

•  In  neutrally  buoyant  systems,  benefits  of  fuel  cell  come  with  minimal  range  cost. 

•  Sensitivity  analysis  illustrates  the  influence  of  design  and  modeling  assumptions. 

•  Predict  up  to  5-  or  8-fold  range  improvement  over  available  battery  technology. 
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A  unique  requirement  of  underwater  vehicles’  power/energy  systems  is  that  they  remain  neutrally 
buoyant  over  the  course  of  a  mission.  Previous  work  published  in  the  Journal  of  Power  Sources  reported 
gross  as  opposed  to  neutrally-buoyant  energy  densities  of  an  integrated  solid  oxide  fuel  cell/Rankine- 
cycle  based  power  system  based  on  the  exothermic  reaction  of  aluminum  with  seawater.  This  paper 
corrects  this  shortcoming  by  presenting  a  model  for  estimating  system  mass  and  using  it  to  update  the 
key  findings  of  the  original  paper  in  the  context  of  the  neutral  buoyancy  requirement.  It  also  presents  an 
expanded  sensitivity  analysis  to  illustrate  the  influence  of  various  design  and  modeling  assumptions. 
While  energy  density  is  very  sensitive  to  turbine  efficiency  (sensitivity  coefficient  in  excess  of  0.60),  it  is 
relatively  insensitive  to  all  other  major  design  parameters  (sensitivity  coefficients  <  0.15)  like  compressor 
efficiency,  inlet  water  temperature,  scaling  methodology,  etc.  The  neutral  buoyancy  requirement  in¬ 
troduces  a  significant  ( ~  15%)  energy  density  penalty  but  overall  the  system  still  appears  to  offer  factors 
of  five  to  eight  improvements  in  energy  density  (i.e.,  vehicle  range/endurance)  over  present  battery- 
based  technologies. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  United  States  Navy  has  identified  certain  key  naval  missions 
like  intelligence,  surveillance,  and  reconnaissance  as  being  best 
performed  by  Unmanned  Undersea  Vehicles  (UUVs)  [1].  Critical  to 
the  success  of  these  missions  are  compact,  quiet,  and  highly  efficient 
power  and  energy  systems  that  enable  long  range  and  endurance  [  1  ]. 
Previous  work  [2-4]  explored  the  viability  of  an  unusual  propulsion 
system  called  the  Hybrid  Aluminum  Combustor  (HAC)  based  on  the 
exothermic  reaction  of  aluminum  with  seawater.  It  used  a  compre¬ 
hensive  thermodynamic  analysis  to  predict  the  overall  efficiency 
and  volume  of  this  system.  These  data  were  used  to  compute  the 
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system’s  effective  volumetric  energy  density  so  that  it  could  be 
compared  to  that  of  other  power  and  energy  systems.  However,  this 
analysis  was  done  without  regard  to  the  system’s  overall  density. 
This  is  important  as  most  underwater  vehicles  must  be  neutrally 
buoyant  [5]  and  making  them  so  usually  requires  adding  empty 
volume  which  lowers  overall  energy  density  -  often  significantly. 
Therefore,  the  objective  of  this  ‘short  communication’  is  to  correct 
previous  estimates  of  the  aluminum-water  propulsion  system’s 
performance  for  the  effects  of  the  neutral  buoyancy  requirement  so 
that  power  and  energy  systems  can  be  compared  on  a  more 
appropriate  baseline.  This  will  be  accomplished  by  updating  earlier 
models  of  a  1000  L  aluminum  combustion-based  power  section 
(appropriate  for  use  in  heavy-weight  to  large  class  UUVs  [1,6])  to 
predict  the  mass  of  individual  components  (in  addition  to  their 
volume)  and  by  introducing  appropriate  ‘dead’  volume  to  ensure 
that  the  overall  system  is  neutrally  buoyant.  This  empty  volume  can 
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Nomenclature 

AHV 

volumetric  energy  content 

At 

elapsed  time,  endurance 

Across 

cross  sectional  area 

£ 

heat  exchanger  effectiveness 

BOP 

balance  of  plant 

Vp 

propulsive  efficiency 

Cd 

drag  coefficient 

Vt 

thermodynamic  efficiency 

Cp 

specific  heat  capacity 

p 

density 

CEA 

chemical  equilibrium  with  applications 

EDv 

volumetric  energy  density 

Subscripts 

HAC 

hybrid  aluminum  combustor 

C 

cold  side  property 

L 

length 

cell 

single  fuel  cell  property 

M 

mass 

comp 

component 

m 

mass  flow  rate 

H 

hot  side  property 

MEA 

membrane  electrode  assembly 

HE 

heat  exchanger 

N 

number/quantity 

inlet 

property  at  component  inlet 

n 

mass  scaling  exponent 

parts 

internal  parts  of  machinery 

NPSS 

numerical  propulsion  system  simulation 

PL 

payload 

Sx 

sensitivity  coefficient  with  respect  to  parameter  X 

reac 

reactant 

SOFC 

solid  oxide  fuel  cell 

reacted 

amount  consumed  by  reaction 

UUV 

unmanned  undersea  vehicles 

ref 

reference  state 

V 

volume 

stack 

fuel  cell  stack 

vc 

cruise  velocity 

sys 

system  property 

W 

power 

also  be  used  to  keep  system  mass  constant  and  maintain  neutral 
buoyancy  by  storing  waste  products  or  water  to  balance  the 
consumed  reactants.  The  effects  of  hydrogen  utilization  (i.eM  the 
percentage  of  H2  generated  by  the  AI-H2O  reaction  that  is  consumed 
by  the  fuel  cell)  and  operating  depth  are  also  re-explored  in  the 
context  of  the  neutral  buoyancy  requirement. 

2.  Propulsion  system  concept 


unless  a  fuel  cell  is  also  present.  Incorporating  a  fuel  cell  requires 
the  following  additional  components:  an  oxidizer  tank  (63%  by 
mass  H2O2),  a  membrane  electrode  assembly  (MEA),  a  condenser 
and  separator  to  trap  and  remove  the  water  produced  in  the  MEA, 
and  a  compressor  to  return  unused  02  to  the  cathode.  A  decom¬ 
position  reactor  is  also  added  to  facilitate  the  breakdown  of  H2O2  to 
H20  and  02  although  it  could  be  possible  to  do  this  directly  on  the 
SOFC’s  cathode. 


Fig.  1  is  a  schematic  illustration  of  the  HAC  system  with  inte¬ 
grated  solid  oxide  fuel  cell  (SOFC)  investigated  in  Ref.  [4].  The  SOFC 
produces  power  in  a  way  that  increases  overall  system  efficiency 
while  enabling  periods  of  all  electric  operation.  It  also  converts  the 
large  volumes  of  waste  H2  produced  by  the  aluminum-seawater 
reaction  to  liquid  water  that  is  much  more  easily  stored.1  The  sys¬ 
tem  operation  is  described  in  detail  elsewhere  [2,4,7]  and  is  only 
briefly  summarized  here:  a  fuel  feeder  produces  a  stream  of 
aluminum  powder  suspended  in  H2  gas  which  reacts  with  steam  to 
form  AI2O3,  more  H2,  steam,  and  heat.  The  flow  passes  through  a 
separator  to  remove  the  aluminum  oxide.  A  portion  (20-25%)  of 
the  flow  is  recirculated  to  sustain  the  reaction  in  the  combustor.  The 
remainder  passes  through  a  turbine  that  extracts  mechanical  po¬ 
wer.  The  flow  passes  to  the  SOFC  where  a  fraction  of  the  hydrogen 
(the  ‘H2  utilization’)  is  consumed  to  produce  electric  power.  A  heat 
exchanger  then  transfers  waste  enthalpy  in  the  exhaust  to  pre¬ 
warm  the  incoming  water.  The  remaining  steam  is  fully 
condensed  and  supplemented  with  seawater  to  make  up  for  what 
was  consumed  by  the  reaction.  A  portion  of  the  H2  gas  is  separated 
from  the  condensed  water,  compressed,  and  reused  for  the 
aluminum  feeder.  Any  remaining  H2  must  be  vented  overboard 


1  Stealth  considerations  often  prohibit  dumping  the  AI2O3  and  H2  waste  pro¬ 
duced  by  the  Al— H20  overboard.  However,  the  thermodynamic  cost  of  compressing 
and  storing  H2  on  board  can  be  large  [7].  The  ready  availability  of  cooling  water 
makes  it  much  simpler  to  achieve  isothermal  compression  —  heat  exchangers  can 
be  smaller,  integrated  into  the  hull,  etc.  —  but  there  is  a  complex  tradeoff  between 
the  efficiency  gain  and  the  volume  displaced  by  these  additional  components. 
Exploring  this  tradeoff  is  beyond  the  scope  of  this  work  and  attention  focuses 
instead  on  identifying  the  limits  of  system  performance  associated  with  adiabatic 
vs.  isothermal  compression. 


3.  Model  formulation 

3.1.  Simulation  environment 

Mathematical  models  are  developed  for  each  component  and 
are  coupled  and  solved  in  an  environment  called  Numerical  Pro¬ 
pulsion  System  Simulation  (NPSS)  [8].  Models  can  range  from 
simple  algebraic  equations  to  links  to  external  thermodynamic 
solvers  like  chemical  equilibrium  analysis  (CEA)  [9  or  computa¬ 
tional  fluid  dynamics  packages.  Additional  information  about  NPSS 
and  the  individual  component  models  is  available  elsewhere  [4,8]. 
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Fig.  1.  Diagram  of  the  HAC-SOFC  hybrid  system  from  Ref.  [4]. 
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The  focus  here  will  be  on  the  modifications  made  to  the  earlier 
model  4]  to  enforce  neutral  buoyancy. 

3.2.  Equations  for  UUV  range  and  endurance 

The  range  and  endurance  of  a  UUV  cruising  at  constant  speed 
(vc)  are  given  by  [7]: 


[7]  to  determine  the  appropriate  length  (I).  The  volume  of  the  heat 
exchanger  also  depends  on  the  required  effectiveness  (r)  and  heat 
capacity  ratio  between  the  streams  (rnHCPjH/mcCp  c). 


m  \  L[c,  (rhHCp;H/thcCp;C)] 
™ref /  t[£ref,  (^HCp,H/^cCp,c)ref] 


(5) 


Range  = 


Vp-EDyV, 


sys 


+  Q  Pseawater^c)  ’  cross) 


(i) 


y-i  1  A  ,  Vp  '  EDy  •  VSys 

Endurance  =  At  =  ^ - — — - -2 -  (2) 

Wpi 

+  (2  Pseawater  vt)-(CDA 

cross) 

In  these  expressions,  Cd  is  the  vehicle  drag  coefficient,  WPL  is  the 
payload  power,  Across  is  the  vehicle  cross-section  area,  Vsys  is  the 
total  volume  of  the  power  and  energy  system,  rjp  is  the  propulsive 
efficiency,  and  EEfy  is  the  ‘effective’  energy  density  of  the  power  and 
energy  system.  The  ‘effective’  energy  density  is  the  total  recover¬ 
able  energy  available  in  the  system  divided  by  the  total  system 
volume: 


EDV  =  VAHv.reac'^  (2 3) 

''sys 

In  this  expression,  AHv,reac  is  the  volumetric  energy  density  of  re¬ 
actants,  Vi-eac/^sys  *s  th e  fraction  of  total  system  volume  allocated  for 
reactant  storage,  and  rjt  is  the  thermodynamic  efficiency  of  the 
power  and  energy  system.  Since  the  objective  is  to  evaluate  the 
aluminum  combustor/SOFC  system  as  a  ‘drop  in’  replacement  for 
existing  power  and  energy  systems,  Vsys  and  are  assumed  to  be 
fixed  here  and  EDy  is  the  performance  parameter  of  interest. 


The  volume  of  the  fuel  cell  is  estimated  based  on  the  required 
number  of  cells  (determined  by  dividing  the  total  required  electric 
power  by  the  design  power  output  of  a  single  cell),  prior  estimates 
of  the  mass  and  volume  of  an  individual  cell  (0.033  kg  and  0.033  L 
for  a  64  cm  cell  [10]),  and  the  assumption  that  the  SOFC  system 
volume  (Vsofc)  including  pipes  and  pumps  is  3  times  the  volume  of 
the  stack  alone  (Vstack)  [10].  Although  this  approach  neglects 
additional  factors  that  could  arise  in  the  case  of  a  very  large  stack 
(which  could  require  proportionally  more  structural  support  and 
less  insulation),  it  provides  a  reasonable  approximation. 

Vjofc  -  (^)  N«„  vcal  =  f%^Vvca,  (6) 

\  Vstack /  stack/  \  Wcell  J 

The  mass  of  a  component  is  assumed  to  scale  with  its  volume  ac¬ 
cording  to  Eq.  (7)  where  the  parameter  n  represents  the  nature  of 
the  scaling.  A  value  of  n  =  1  corresponds  to  a  linear  scaling  of 
component  mass  with  volume.  This  is  valid  for  thin  walled  cylin¬ 
drical  pressure  vessels  where,  for  a  fixed  operating  pressure,  wall 
thickness  is  proportional  to  the  radius  [11  ]  and  therefore  wall  mass 
is  proportional  to  radius  squared  times  length  (i.e.,  volume).  A  value 
of  n  =  2/3  corresponds  to  a  situation  where  mass  scales  roughly 
with  surface  area.  This  is  valid  for  a  vessel  with  constant  wall 
thickness.  The  influence  of  the  choice  of  scaling  exponent  will  be 
examined  by  comparing  results  obtained  using  each  assumption. 


3.3.  Volume  and  mass  scaling 

The  volume  of  each  component  is  assumed  to  be  directly  pro¬ 
portional  to  the  flow  passing  through  it  so  that  components  get 
larger  or  smaller  with  the  design  power  requirement  of  the  system. 
Each  component  is  assumed  to  be  cylindrical,  flow  velocities 
through  them  are  assumed  to  be  uniform,  and  residence  time  is 
assumed  to  be  constant.  These  assumptions  lead  to  Eq.  (4)  which 
describes  how  the  volume  of  a  component  varies  with  respect  to  a 
reference  size  associated  with  the  baseline  system.  Reference  states 
for  various  components  consistent  with  a  15  kW  system  are  listed 
in  Table  1. 

Vcomp  =  (4) 

Heat  exchangers  are  assumed  to  be  simple  straight  tube,  counter- 
current  designs.  An  iterative  finite  differencing  scheme  is  applied 


Table  1 

HAC  system  component  reference  states. 


Component 

mref(kgs  ’) 

l^comp.ref  (L) 

‘Pack  factor’ 

Mcomp.ref  (kg) 

Fuel  feeder 

0.0147 

6.25 

0.25 

5.0 

Combustor 

0.0835 

50 

0.05 

14.2 

Oxide  separator 

0.0835 

25 

0.15 

18.1 

Steam  compressor 

0.0147 

6.25 

0.35 

10.1 

Steam  turbine 

0.0571 

25 

0.35 

40.3 

Heat  exchanger 

0.1159 

62.5 

0.30 

63.2 

Water  pumps 

0.0587 

25 

0.35 

23.8 

Condenser 

0.0571 

20 

0.20 

10.8 

H2  separator 

0.0571 

25 

0.15 

10.2 

H2  compressors 

0.0014 

5 

0.35 

4.8 

Mcomp  _  /  Vcomp 
^comp,ref  yVCOmp,ref 

The  mass  of  the  reference  configuration  (Table  1 )  was  estimated 
using  the  following  very  basic  sizing  techniques.  The  component 
wall  thicknesses  were  selected  to  withstand  pressure  stresses  using 
standard  techniques  for  cylindrical  pressure  vessel  design  [11]  and 
a  factor  of  safety  of  2.  The  remaining  mass  is  determined  by  making 
some  basic  assumptions  about  the  internal  design  of  turbines, 
compressors,  heat  exchangers,  etc.  in  order  to  estimate  a  ‘packing 
factor’  equal  to  the  fraction  of  interior  volume  filled  by  structural 
elements  or  mechanical  parts  (=Vparts/Vc omp)-  A  ‘packing  factor’  of  1 
corresponds  to  solid  block  whereas  0  corresponds  to  a  hollow  shell. 
The  ‘packing  factors’  used  in  this  work  are  listed  in  Table  1.  The 
largest  factors  are  for  the  pumps,  turbine,  and  compressors  which 
have  lots  of  moving  parts.  The  smallest  factor  is  for  the  combustor 
which  is  a  nearly  hollow  chamber.  Most  components  are  assumed 
to  be  constructed  from  aluminum  but  titanium  alloys  are  assumed 
for  high  temperature  components  like  the  combustor  and  turbine. 

3.4.  Energy  density  calculations 

The  system  volume  is  the  sum  of  the  HAC,  SOFC,  reactants,  and 
any  empty  volume  needed  to  achieve  neutral  buoyancy.  The  mass 
available  for  reactant  (fuel)  storage  is  found  by  requiring  that  the 
total  mass  of  the  system  equals  the  volume  of  the  system  times  the 
density  of  water.  This  leads  directly  to  an  expression  for  the  volume 


2  Note  that  the  47%  number  depends  on  the  performance  of  the  aluminum  seeder 

and  could  be  different  for  different  systems. 
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of  fuel  required,  Eq.  (8).  Rewriting  the  product  of  efficiency  and 
reaction  enthalpy  from  Eq.  (3)  as  the  ratio  of  system  power  to 
reactant  flow  rate  gives  Eq.  (9)  for  the  system  energy  density  in 
terms  of  the  quantities  in  the  model. 

u  _  P water ^sys  ~  MHac  ~  MS0FC 


EDy  =  (Whac  +  Wsofc)|^^ 

4.  Results  and  discussion 

4.1.  Operating  cases  considered 

The  H2  utilization  rate  (=  mH2  reacted/mH2  inlet)  is  set  to  0%  for 
the  first  model  run  in  order  to  establish  a  baseline  corresponding  to 
a  HAC  system  operating  without  a  fuel  cell.  The  full  range  of 
available  H2  utilizations  is  examined  from  0%  up  to  53%  where  all  of 
the  ‘excess’  H2  is  consumed.  The  remaining  47%  of  the  H2  produced 
by  the  A1-H20  reaction  is  needed  to  operate  the  aluminum  feeder. 
In  this  analysis,  any  ‘excess’  H2  is  assumed  to  be  vented  overboard 
and  the  work  required  to  compress  it  to  the  ambient  pressure  (so 
that  it  can  be  dumped)  is  accounted  for  when  computing  the  net 
power  output  of  the  system.  Depths  of  3  m  (131.7  kPa),  30  m 
(407.5  kPa),  and  150  m  (1633.4  kPa)  are  investigated  to  cover  the 
operating  range  defined  as  necessary  by  the  Navy  for  UUVs  in  this 
class  [6,12]. 

Calculations  for  each  case  begin  at  0%  H2  utilization  and  15  kW 
net  power  output.  This  power  level  is  chosen  in  order  to  be 
consistent  with  a  1000  L  power  section  operating  at  10-20  W  L  1 
power  density  [12].  The  HAC  system  is  scaled  as  required  by  varying 
the  fuel  mass  flow  (which  in  turn  causes  other  flows  in  the  system 
to  vary).  The  efficiencies  of  all  components  are  assumed  to  remain 
constant  over  the  duration  of  the  mission. 

4.2.  Effects  of  H2  utilization  and  depth 

Fig.  2  shows  energy  density  as  a  function  of  fuel  cell  H2  utili¬ 
zation  assuming  adiabatic  (lower)  and  isothermal  (upper  curves) 
H2  compression  at  70%  adiabatic  and  isothermal  efficiency, 


(8) 

(9) 


respectively.  Recall  that  adiabatic  and  isothermal  compressions  are 
assumed  to  identify  the  lower  and  upper  bounds  of  system  per¬ 
formance  respectively.  The  thick  curves  correspond  to  the 
assumption  that  HAC  mass  scales  linearly  with  volume,  and  the 
thin  curves  just  below  them  correspond  to  the  assumption  that 
mass  is  proportional  to  surface  area.  The  linear  scaling  leads  to 
higher  energy  densities  than  predicted  by  area  proportional  scaling. 
However,  the  difference  is  relatively  small  (1%— 3.5%)  so  the  linear 
scaling  is  used  for  the  rest  of  the  calculations  reported  here. 

With  adiabatic  compression  at  3  m,  the  system  delivers  a 
neutrally  buoyant  energy  density  of  277  W  h  L_1  at  0%  H2  utilization 
(i.e.,  the  baseline  HAC  system)  and  509  W  h  L_1  at  full  (53%)  utili¬ 
zation.  The  neutrally  buoyant  energy  density  peaks  at  532  W  h  L-1 
and  35%  H2  utilization  because  of  a  competition  between  system 
efficiency  (which  is  improved  by  incorporating  the  fuel  cell)  and  the 
displacement  of  fuel  by  the  fuel  cell  system.  With  isothermal 
compression  at  3  m,  the  system  delivers  a  neutrally  buoyant  energy 
density  of  714  W  h  L-1  (the  maximum  observed)  at  0%  H2  utiliza¬ 
tion  and  606  W  h  L-1  (15%  below  peak)  at  full  utilization.  No  peak  is 
observed  because  of  the  improved  H2  compression  efficiency. 
Switching  to  isothermal  compression  more  than  doubles  the  en¬ 
ergy  density  of  the  system  at  zero  H2  utilization.  Isothermal 
compression  improves  energy  density  by  19%  at  full  utilization. 
While  these  results  suggest  that  increasing  system  complexity  by 
using  cooled  compressors  could  be  very  worthwhile,  this  only 
represents  an  upper  bound  on  what  is  possible  thermodynamically 
because  no  accounting  has  been  made  for  the  additional  volume 
required  and  the  performance  of  a  real  70%  efficient  system  would 
be  less.  The  effects  of  increasing  vehicle  depth  are  also  clear  in 
Fig.  2:  As  depth  is  increased,  more  work  is  required  to  raise  the 
hydrogen  to  the  ambient  pressure  which  reduces  system  energy 
density.  The  curves  converge  at  full  H2  utilization  where  there  is  no 
need  to  vent  excess  H2  and  the  system  becomes  depth  independent. 
However,  if  stealth  is  required  (i.e.,  venting  H2  overboard  is  pro¬ 
hibited)  then  the  best  performance  always  occurs  at  100%  H2 
utilization. 

Fig.  3  compares  the  neutrally  buoyant  results  to  those  previously 
reported  [4]  for  negatively  buoyant  (denser  than  water)  systems. 
Only  results  at  3  m  depth  are  shown  (recall  that  surface  area  pro¬ 
portional  mass  scaling  results  will  be  1-3.5%  lower).  The  compar¬ 
ison  shows  that  the  penalty  associated  with  enforcing  neutral 
buoyancy  is  quite  large:  energy  density  is  reduced  by  58%  at  zero  H2 


Percent  H2  Utilization  Percent  H2  Utilization 


Fig.  2.  System  energy  density  vs.  hydrogen  utilization  under  neutrally  buoyant  con-  Fig.  3.  System  energy  density  vs.  hydrogen  utilization,  comparison  of  neutral  and 
ditions  at  various  depths.  negative  buoyancy. 
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Fig.  4.  Component  volume  vs.  hydrogen  utilization  for  neutrally  buoyant  system  at 
3  m  depth.  Isothermal  compression  is  assumed. 


utilization  and  by  29%  at  full  utilization  for  either  the  adiabatic  or 
isothermal  compression  cases.  The  penalty  is  smaller  at  high  uti¬ 
lization  because  the  fuel  cell  and  H2O2  are  less  dense  than  the  HAC 
and  aluminum  they  displace.  This  is  illustrated  in  Fig.  4  which 
shows  the  volume  distribution  in  the  1000  L  power  and  energy 
section  as  a  function  of  H2  utilization  for  the  neutrally  buoyant 
system  at  3  m  depth  using  isothermal  compression.  The  HAC,  SOFC, 
and  reactant  storage  volumes  are  shown  along  with  the  empty 
volume  necessary  to  achieve  neutral  buoyancy.  The  small  cross- 
hatched  volume  illustrates  the  difference  between  the  predictions 
of  the  two  mass  scaling  assumptions.  Assuming  HAC  mass  pro¬ 
portional  to  volume  predicts  a  slightly  lighter  system  which  allows 
for  more  reactant  storage  and  thus  higher  overall  system  energy 
density.  The  figure  shows  that  the  volume  that  must  be  left  empty 
at  zero  utilization  is  approximately  twice  that  required  at  full  uti¬ 
lization  and  so  the  penalty  associated  with  adding  the  fuel  cell  is 
reduced.  The  net  result  is  that  for  isothermal  compression  at  a 
depth  of  3  m,  adding  the  fuel  cell  only  drops  the  neutrally  buoyant 
energy  density  by  15%  compared  to  50%  in  the  denser  than  water 
system.  Therefore,  the  benefits  of  adding  a  fuel  cell  (eliminating  H2 
venting,  depth  independence,  enabling  periods  of  all  electric 
operation)  come  at  much  smaller  cost  when  the  effects  of  main¬ 
taining  neutral  buoyancy  are  factored  in. 


4.3.  Sensitivity  analysis 

Since  so  many  assumptions  must  be  made  in  analyses  of  this 
type,  it  is  useful  to  quantify  the  impact  of  the  impact  of  these  as¬ 
sumptions  on  the  predicted  energy  density.  This  is  accomplished  by 
computing  ‘sensitivity  coefficients’  (Sx)  for  each  parameter  X.  It  is 
defined  as  the  fractional  change  in  system  energy  density  divided 
by  the  fractional  change  in  the  parameter  as  shown  in  Eq.  (10) 
where  ‘re/  refers  to  the  baseline  assumptions  and  results  and 
prime  values  denote  the  values  of  the  perturbed  parameter  and  the 
resulting  energy  density. 


Sx  = 


(ED'V  EDy  ref^  / ED\], ref 


X'  -X, 


ref 


/X, 


(10) 


ref 


The  results  of  the  sensitivity  analysis  are  shown  in  Fig.  5.  All 
changes  are  relative  to  the  neutrally  buoyant  HAC  system  (HAC 
mass  aV)  with  an  isothermal  compression  system  operating  at  3  m 
depth  and  full  utilization.  The  various  parameters  investigated  are 
grouped  into  four  categories:  component  performance  parameters 
(turbine  efficiency,  compressor  efficiencies,  fuel  seed  ratio),  fuel  cell 
material  properties  (electrode  porosity  and  tortuosity,  thickness  of 
electrodes  and  electrolytes),  system  operating  parameters  (tem¬ 
peratures,  pressures,  current  density,  power  demand),  and  scaling 
assumptions  (reference  masses  and  volumes). 

The  figure  shows  that  the  turbine  efficiency  has  by  far  the 
largest  impact  on  the  overall  system  energy  density  and  thus  is  the 
most  important  parameter  to  know  well.  The  value  used  here  (0.60) 
is  well  below  typical  values  for  large,  state  of  the  art  steam  turbines. 
The  value  is  chosen  to  be  a  conservative  estimate  for  smaller  sized 
machines  in  this  type  of  application  [13].  Another  particularly 
important  parameter  is  the  H2  compressor  efficiency.  The  value 
used  here  (0.70)  is  a  conservative  estimate  for  typical  hydrogen 
compressors  [14]. 

Several  system  operating  parameters  also  have  important  in¬ 
fluences  on  cycle  performance,  but  unfortunately  the  operator  has 
only  limited  ability  to  control  many  of  them.  The  incoming  water 
temperature  is  set  by  environmental  conditions,  turbine  pressure 
ratio  is  limited  by  system  pressures,  and  power  demand  is  set  by 
mission  requirements.  Note  that  the  assumed  inlet  water  temper¬ 
ature  (21  °C)  is  relatively  high  (allowing  for  seasonal  and  regional 
differences),  and  operating  at  temperatures  that  are  more  repre¬ 
sentative  of  ocean  averages  would  substantially  improve  perfor¬ 
mance.  The  results  also  show  that  it  would  be  beneficial  to  raise 
overall  system  pressures,  but  this  will  be  limited  by  the  need  to 
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increase  wall  thicknesses,  etc.  to  withstand  higher  pressures.  The 
additional  mass  that  this  would  require  is  not  accounted  for  in  the 
current  model.  The  properties  of  the  materials  used  to  construct  the 
fuel  cell  also  have  significant  impacts  on  energy  density,  but  their 
importance  is  at  least  an  order  of  magnitude  smaller  than  the 
turbine  efficiency.  For  example,  the  electrodes’  porosity  and  tor¬ 
tuosity  control  diffusion  and  thus  the  diffusion  overpotential. 
Anode  thickness  has  more  impact  than  cathode  or  electrolyte 
thickness  because  of  the  importance  of  anode  diffusion  in  this 
particular  SOFC  arrangement.  The  values  of  these  parameters  are 
consistent  with  those  used  in  other  SOFC  modeling  studies  [15-18]. 
Sensitivities  to  choices  of  SOFC  reference  mass,  volume,  and  bal¬ 
ance  of  plant  (BOP)  volume  are  also  approximately  one  order  of 
magnitude  smaller  than  the  sensitivity  to  turbine  efficiency.  The 
volume  of  the  SOFC  balance  of  plant  components  is  related  to  the 
V^soFe/^stack  parameter  [10]  discussed  in  Section  3.3  (i.e.,  Vbop / 
Vstack  =  V^soFe/^stack  -  1).  The  choice  of  SOFC  reference  mass  and 
volumes  have  larger  impacts  than  the  choice  of  HAC  reference 
volume  and  mass  because  at  full  utilization  the  SOFC  is  the  larger 
system. 

Boundaries  on  the  overall  (cumulative)  uncertainty  in  the  pre¬ 
dicted  energy  density  were  estimated  by  perturbing  all  of  the 
design  parameters  by  10%  in  a  coordinated  manner  so  that  they 
either  raise  or  lower  the  overall  system  energy  density.  The  results 
show  that  a  ±10%  shift  in  these  parameters  translates  to  a  ±15% 
shift  in  predicted  energy  density.  Since  it  is  unlikely  that  the  values 
of  all  of  the  design  parameters  would  be  off  by  10%  and  that  the 
deviations  would  all  bias  the  energy  density  in  the  same  direction, 
15%  most  likely  represents  an  upper  bound  on  the  uncertainty  of 
the  energy  density  estimates. 

4.4.  Overall  performance  comparison 

Fig.  6  compares  predicted  energy  densities  for  FIAC  and  HAC- 
SOFC  systems  to  those  of  other  UUV  power  systems  based  on 
SOFCs  [10],  lithium-ion  [19-21],  lead-acid  [21],  and  silver-zinc  [21] 
battery  systems.  Results  for  both  the  present  (neutrally  buoyant, 
both  scaling  assumptions)  and  past  (negatively  buoyant  [4])  studies 
are  presented  for  comparison.  The  ‘error’  bars  indicate  the 
approximate  range  of  performance  that  one  could  expect  based  on 
the  reports  in  the  literature.  For  the  HAC  systems,  the  solid  portion 
of  the  bars  shows  the  performance  with  adiabatic  compression.  The 


very  top  of  the  cross-hatched  bars  shows  the  performance  that  is 
possible  with  an  isothermal  compressor  that  adds  negligible  mass 
and  volume  compared  to  an  adiabatic  compressor.  The  intermedi¬ 
ate  lines  in  the  cross-hatched  portion  show  the  peak  performance 
for  isothermal  compression  systems  of  various  masses  and 
volumes. 

The  figure  shows  that  six-  or  seven-fold  improvements  in 
neutrally  buoyant  system  energy  density  are  possible  over  existing 
Li-ion  battery  technology  with  compact  isothermal  compression. 
While  the  neutral  buoyancy  requirement  results  in  a  significant 
drop  in  energy  density,  battery-based  systems  suffer  similar 
penalties.  The  pure  SOFC  system  fares  better  under  the  neutrally 
buoyant  restriction  because  the  system  defined  in  Ref.  [10]  is 
already  nearly  neutrally  buoyant.  Note  that  it  only  slightly 
underperforms  the  peak  performance  of  the  HAC  and  HAC-SOFC 
systems. 

The  neutral  buoyancy  constraint  also  greatly  reduces  perfor¬ 
mance  differences  between  the  various  HAC  cases  compared  to  the 
earlier  results  in  the  heavier  than  water  systems.  This  is  an 
important  finding  because  the  earlier  results  suggested  that  elim¬ 
inating  gas  venting  and  achieving  depth  insensitivity  came  at  very 
large  cost  in  terms  of  reduced  performance.  The  present  results 
show  that  in  neutrally  buoyant  systems,  much  less  performance 
must  be  sacrificed  in  order  to  achieve  the  important  benefits 
offered  by  the  HAC-SOFC  hybrid  system. 

5.  Conclusions 

An  earlier  numerical  model  of  an  underwater  power  system  has 
been  expanded  to  account  for  the  fact  that  most  underwater  power 
and  energy  systems  must  be  neutrally  buoyant.  While  this  imposes 
a  substantial  energy  density  penalty,  other  systems  suffer  similar 
penalties  and  the  hybrid  aluminum  combustor/solid  oxide  fuel 
system  fares  proportionally  better.  Another  key  finding  is  that  the 
neutrally  buoyant  systems  show  far  less  variation  between  the 
straight  aluminum  burning  system  and  the  fuel  cell  hybrid  system. 
This  means  that  the  benefits  of  the  fuel  cell  (elimination  of  H2 
venting,  depth  independence,  and  enabling  periods  of  all  electric 
operation)  are  achieved  with  minimal  performance  penalty.  Finally, 
a  sensitivity  analysis  shows  that  the  performance  of  the  turbine  has 
by  far  the  largest  impact  on  the  overall  energy  density  achievable 
by  the  system.  The  baseline  model  assumptions  suggest  that  energy 
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Fig.  6.  Comparison  of  the  energy  densities  of  various  underwater  power  and  energy  technologies. 
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density  improvements  of  factors  of  6-7  over  current  technology 
(batteries)  are  possible.  This  range  broadens  to  factors  of  5-8  im¬ 
provements  in  overall  system  energy  density  when  uncertainties  in 
the  model’s  assumptions  and  the  system’s  sensitivity  to  them  are 
accounted  for.  Therefore,  the  aluminum  combustion  system  ap¬ 
pears  to  be  an  extremely  promising  power/energy  source  for  un¬ 
manned  underwater  vehicles. 
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